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l" ■ In the Horndeski's most general scalar-tensor theories the equations of scalar density perturba- 

I ' tions are derived in the presence of non-relativistic matter minimally coupled to gravity. Under a 

, quasi-static approximation on sub-horizon scales we obtain the effective gravitational coupling G e ff 

■ associated with the growth rate of matter perturbations as well as the effective gravitational poten- 
tial t&off relevant to the deviation of light rays. We then apply our formulas to a number of modified 
gravitational models of dark energy-such as those based on f(R) theories, Brans-Dicke theories, ki- 
netic gravity braidings, covariant Galileons, and field derivative couplings with the Einstein tensor. 

' Our results are useful to test the large-distance modification of gravity from the future high-precision 

\ observations of large-scale structure, weak lensing, and cosmic microwave background. 

<n : 

i— i' I. INTRODUCTION 

O . 

cr ... , 

The late-time cosmic acceleration has been supported by several independent observations-such as supernovae la 
tjjTj' [H> Cosmic Microwave Background (CMB) 0, and baryon acoustic oscillations Q. The simplest candidate for dark 
i— i, energy is the cosmological constant, but the typical scale of the vacuum energy is vastly larger than the observed 
energy scale of dark energy Instead many alternative models have been proposed to identify the origin of dark 
<^j; energy 0,0. 

A minimally coupled scalar field with a potential V(</>)-quintessence [jj-can account for the cosmic acceleration 
today, provided that the potential is sufficiently flat with a small effective mass « 10~ 33 eV. The k-essence 
, where the Lagrangian includes a nonlinear term of the field kinetic energy, can be responsible for dark energy even in 
the absence of a field potential. Quintessence and k-essence can be distinguished from the cosmological constant in 
qq , that their equations of state vary in time while the latter does not. 

■ There is another class of dark energy models based on the large-distance modification of gravity-such as (i) f{R) 
t— I \ theories , (ii) Brans-Dicke theories (hi) Dvali-Gabadadze-Porrati (DGP) braneworld [ll|, and (iv) Galileon 
t-H . gravity [12 1. In the local region where the average density is much larger than the cosmological one, these models 

need to recover the General Relativistic behavior for consistency with solar-system experiments [l3[. The models 
• , based on the theories (i) and (ii) can be made to be compatible with local gravity constraints under the chameleon 
mechanism [l4| , as long as the scalar degree of freedom has a large effective mass in the region of high density [ll| . 
For the models based on the theories (iii) and (iv) the nonlinear field self-interaction can allow the recovery of the 
General Relativistic behavior in the local region [To ] through the Vainshtein mechanism [l7j . 

For the dark energy models mentioned above the field equations of motion are kept up to second order. This is 
desirable to avoid the appearance of the Ostrogradski's instability pj|] associated with the derivatives higher than 
the second order. In 1974 Horndeski pj| derived the most general single-field Lagrangian for scalar-tensor theories 
with second-order equations of motion. Recently this issue was revisited by Deffayet et al. [20] in connection to a 
covariant Galileon field. The most general scalar-tensor theories with the second-order equations can be expressed by 
the sum of the Lagrangians (0-([5]) below. In fact one can show that this Lagrangian is equivalent to that derived by 
Horndeski [2l[ (see also Ref. |22j). 

The most general scalar-tensor theories not only include quintessence and k-essence but also accommodate f(R) 
theories, Brans-Dicke theories, and Galileon gravity. Moreover, as shown in Ref. [2l[ ,several different choices of the 
functions A', Gi (i = 3,4,5) give rise to the (modified) DGP model in 4 dimensio ns 12311 . the field coupling with the 
Gauss-Bonnet term [24( , the field-derivative coupling with the Einstein tensor [25l . |26|. and so on. 

In this paper we shall derive the equations of linear density perturbations for the most general scalar-tensor theories 
with non-relativistic matter taken into account. In the presence of the terms Ci (i = 3, 4, 5) the effective gravitational 
coupling G ff is subject to change compared to that in General Relativity. This leads to the modified growth rate of 
matter density perturbations S rn as well as the modified evolution of the effective gravitational potential $ c ff associated 
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with the deviation of light rays. Similar analysis has been carried out in specific scalar-tensor theories (6l [27l - l3l| . f(R) 
theories (3ll . |32T | , kinetic gravity braidings with the term C 3 f33T - [37| , and covariant Galileon [38| . Our analysis in this 
paper covers those theories as specific cases. Such general analysis will be useful to discriminate between modified 
gravitational models from the observations of large-scale structure, weak lensing, and CMB (39j . 

This paper is organized as follows. In Sec. [II] we derive the background equations of motion on the flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) background for the action ([6]) below. In Sec. Mil we obtain the full scalar 
perturbation equations of motion for the metric (|21|) . In Sec. IIVI the effective gravitational coupling G e e as well as 
the effective gravitational potential <E» c ff are derived under the quasi-static approximation on sub-horizon scales. In 
Sec. [V] we apply our formulas of sub-horizon perturbations to a number of modified gravitational models of dark 
energy. Sec. IVIl is devoted to conclusions. 



II. THE MOST GENERAL SCALAR-TENSOR THEORIES AND THE BACKGROUND EQUATIONS 

OF MOTION 

The most general 4-dimensional scalar-tensor theories keeping the field equations of motion at second order are 
described by the Lagrangian [20| 



£ = E £ " (!) 



i=2 



where 



£ 2 = K{4>,X), (2) 
C 3 = -G 3 ((t>,X)n<t>, (3) 
U = G 4 (0,X)i? + G4,x[(n^) 2 -(V A1 V^)(V' I V»], (4) 

C 5 = G 6 (4>, X) (VV» - i G 5 , x [(D0) 3 - 3(D0) (V^V^) (V^V» + 2(V^V Q 0) (V Q Vtf) (V^V M 0)] . (5) 

Here K and Gi (i = 3,4,5) are functions in terms of a scalar field and its kinetic energy X = —d fJ, <fid^(f>/2 with 
the partial derivatives G^x = dGi/dX, R is the Ricci scalar, and G M „ is the Einstein tensor. The above Lagrangian 
was first discovered by Horndeski in a different form [ljj. In fact the Lagrangian JTJ) is equivalent to that derived by 
Horndeski [2l|. 

We are interested in the late-time cosmology in which the field <f> is responsible for dark energy. In addition we 
take into account a barotropic perfect fluid with the equation of state w = P m /Pm, where P m is the pressure and p m 
is the energy density respectively. In the following we focus on non-relativistic matter (w = 0) minimally coupled to 
the field 0. 1 The total action we are going to study is then given by 

S = J d 4 x^(C + C m ) , (6) 

where g is a determinant of the metric g^ u , and C m is the Lagrangian of non-relativistic matter. 

Let us consider a flat FLRW background with the metric ds 2 = —N 2 (t)dt 2 + a 2 (t)dx 2 . Variations with respect to 
the lapse N(t) and the scale factor a(t) give rise to the following equations of motion respectively 



£ = £i = ~Pm , (7) 

»=2 
5 



i=2 



1 If matter is non-minimally coupled to <f> through the coupling to the metric g M „ = A{(j>)gfiu rather than g^v, one rewrites the action JHJ 
in terms of (f> and <? M „ rather than and g^v Then, the Lagrangian is still of the form (0-J5}, because the change of the variable 
9tu> 9^v — A(4>)g^ v does not generate higher derivative terms in the field equations. For this reason, it is sufficient to consider matter 
minimally coupled to <f>. 
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and 



2 HX + HX + 3H 2 X 



where 

£2 = 2XK, X -K, (9) 

S 3 = BX4>HG3,x-2XG3j, (10) 

£ 4 = -6H 2 G 4 + 24H 2 X(G4,x+XG 4 ,xx)-12HX<i>G 4 , 4> x-QH<i>G4,<i > , (11) 

S 5 = 2H 3 X<P {5G 5 , X + 2XG 5 ,xx) - 6H 2 X (3G 5 , + 2XG^ X ) , (12) 

V 2 = K, (13) 

V3 = -2x(g 3 , 4 , + ^G 3 , x ) , (14) 
Vi = 2 (SH 2 + 2Hj G 4 - 12H 2 XGa, x - AHXd.x - SHXG^x - 8HXXG 4 ,xx + 2 + 2H<p) G 4> 

+4XG 4 , 00 + 4X U-2H^j G^x, (15) 
V 5 = -2X (2H 3 4> + 2HH4> + 3H 2 ^j G 5 , x - 4H 2 X 2 4> G 5 , X x + 4HX (x - HX^j G^x 

Gs,t + AHXj>Gw (16) 

Here a dot represents a derivative with respect to t and H = a/ a is the Hubble parameter. Varying the action (JB]) 
with respect to 4>(t), it follows that 

where 

J = 0X x + 6 J ffXG 3 ,x-2(/)G3,0 + 6 J ff 2 0(G4,x + 2XG4 : xx)-12i7XG4,^Y 

+2H 3 X (3G 5 , X + 2XG 5 ,xx) - 6H 2 (G M + XG^x) , (18) 
P = - 2X (g zm + 4>G 3 ,4x) + 6 (2H 2 + ff) G 4 , + 6H (x + 2HX) G^ x 

-6H 2 XG 5 ^ + 2H 3 X<P G^x ■ (19) 
Non-relativistic matter obeys the continuity equation 

p m + 3H Pm = . (20) 

Equations (O, ©, (flTf . and ([20]) are not independent because of the Bianchi identities. In fact the field equation 
(fT7)) can be derived by using Eqs. 0, flSJ), and (|2U|) . 

III. LINEAR PERTURBATION EQUATIONS 

In this section we derive the linear perturbation equations for the theories give by the action (J5]). Let us consider 
the following perturbed metric about the flat FLRW background [4(| 

ds 2 = -(1 + 2*) dt 2 - 2d i xdtdx i + a 2 (t)(l + 2$)S ij dx i dx j , (21) 

where 'J', $, and x are the scalar metric perturbations. In this expression we have chosen a spatial gauge such that 
gij is diagonal, which fixes the spatial part of a vector associated with a scalar gauge transformation. The temporal 
part of the gauge-transformation vector is not fixed for the moment. 

We perturb the scalar field as 4>(t)+S(j)(t, x), and the matter fields as well, in terms of the matter density perturbation 
dp m , and the scalar part of the fluid velocity v. We define the density contrast of matter as S = Sp m /p m . In order to 
write the perturbation equations in a compact form, we introduce the following quantities (2l| . 



F T = 2 
Qt = 2 



G 4 -x(0G5.x + G 5 , )] , (22) 
G 4 - 2XG 4 , X - X (Hi G 5 , x - G M ) ] , (23) 
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and 
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£ = 



1 fig 
~6dH 

-<pXG 3 , x + 2HG 4 - 8HXG^ X - 8HX 2 G 4 , X x + 0G 4 . + 2X<j>G^x 
-H 2 <j> {hXG^x + 2X 2 G 5 , XX ) + 2HX (3G 5 ^ + 2XG^ X ) , 

X d£ + l H 9E 

XK X + 2X 2 K^xx + 12H<j)XG 3 .x + 6H<j)X 2 G 3 ,xx — 2XG 3 — 2X 2 G 



3,4>X 



-6H 2 G 4 + 6 



H 2 (7XG 4tX + WX 2 G 4 ,xx + AX 3 G 4 ^xxx) — H<j) (G4 t< p + bXGi^x + 2X 2 G 4t 



4>xx j 



(24) 



+30H 3 <])XG 5 ,x + 26H 3 <j)X 2 G 5 ,xx + 4H 3 (j)X 3 G5 t xxx — QH 2 X (66*5^ + 9XGf,^x + 2X 2 Gf,^xx) , (25) 

where we used the relation Xdx<t> = 4>/2. The functions Tt and Gt appear in the quadratic action for the cosmological 
tensor perturbations (2ll l4ll , [43 | . In order to avoid ghost and Laplacian instabilities in the tensor sector we require 
the conditions Tt > and Qt > 0. 

We expand the action (|6]) up to second order in perturbations and vary the second-order action with respect to each 
perturbed variable such as ty. Following the procedure explained in Ref. (38j . the perturbation equations in Fourier 
space are given by 



Eq? 

Escj, 

E v 
E d 



k 2 k 2 ( k 2 

A x § + A 2 5(j) - p m v + A 3 —<Z> + Ai* + A 5 — x + A 6 — - fx 

a* a z \ a z 

k 2 ( k 2 
- B^> + B 6 ^$ + B 7 ^r + 3v 



B 4 <P + B 2 5(j) + B 3 $ + B 4 
k 2 



k 2 k 2 
B R —+B 9 )m + B w —x + B 11 —x + 3p m v = 0, 



Ci* + C 2 S(j) + C3* + C 4 8<j) + p m v = , 



+ D 3 <S> + D 4 



D 



D 4 $ + D 2 t 
k 2 

v - * = , 



k 2 k 2 

(5 + 3$ + T w 5-x = 0, 

a a 



D 5 V + D 6 — X 



1 ^ 



9 ~ 



-M 2 )5<j>+(D 



- p m S = , 



k 2 

D n ) * + D 12 — x = 0, 



(26) 

(27) 
(28) 

(29) 
(30) 
(31) 



where k is a comoving wavenumber, and 



Ai=69, A 2 = -2(E + 3H&)/<t>, 
A 6 = 2(9 - Hg T )/<p, » = 



a 3 = 2g Tl 



A 4 = 2Z + p m , A 5 = -26, 



(32) 



B 1= 6£ T , B 2 = 6(Q-Hg T )/j>, B 3 = 6(Sr + 3Hg T ), 



B 4 = 3 
Be = 2Tt 



[AH<j) - AH<f) - 6H (f> ) Q T - 2H<j) Q T 



6Hcf)) 9 + 206-/^ 



B, 



B 7 = 2 



g T + H(g T -T T ) /<t>, B s = 2g Tl b 9 = -6(e + 3HQ), 



B w = -2g T , B n = -2(g T + Hg T ), v = v A 



-66, 



(33) 



C 4 = 2g T , C 2 = 2{Q-Hg T )/4>, C 3 = -2Q, C 4 = 2(H ifi — H (j))g T — 24> Q — p m (j) /(f 2 , (34) 
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£>i = 6(9 - HG T )/<j>, D 2 = 2(3H 2 G T - 6H0 - £)/0 2 , 



D 3 



2H(g T + 3>HQ T ) - 2(9 + 3HG) - p m J /0, 
3H{(3H 2 + 2H)4> - 2H^}g T + 3H 2 j>g T + 6{2i?0 - (3H 2 + H)<j>}Q - 6#09 + (20 - 3#0)£ - <pt] /<fi 



D 5 = 2(S + 3HQ)/j>, D 6 = -2(e-Hg T )/j>, D 7 = 2 Qt + H {Q T - T T ) /<j>, 



D w 



= 3 6(H<j> - H4>)Q - 20T, + 3Hp m 4> - p4> 2 /4> 2 , 



2H 2 T T - 4H(g' T + Hg T ) + 2(9 + HQ) 



I4>\ 



D 1 

M 



2(Q- Hg T )/<j>, D u = 6{(3H 2 + H)<p- Hcf>}Q + 6H<]>Q + 2(3H<p- 0)E + 2</>E - n<j) 2 /<£ 
2H(g T + Hg T ) - 2(9 + HQ) - p m ] /4> , 
[ft + 3H(p + v)] /4> 

-K M + {4> + 3H<P)K^ X + 2XK MX + 2X4>K t<t>X x 
+ [6 H(Gz^xxX + Gz^x)4> — 2 G^^xX — 2 G^^\(j) + 6 H (G^^xX 



G s 



G 



3,000 



\X 



*4,<t>X 



6 G 3:tj>x XH + 2(9 H 2 G 3 ,<j>X 

[6 7J 2 (4 Gi )( f,xxxX 2 + 8 G^^xxX + G^x) — 6 H (2 G^^xxX + 3 G^^x, 
[12 H(G±^x + 2 G4 lP xxX)H + 6H(6H 2 G A ,<f>xxX — 2 Gi^^xX + 3H G4 
12 H 2 (2 Gi^xxX 2 — 3 Gi^xX — G^cf,) — 6 (2 G^^xX + G±^$) H 
[2 H 3 (2G5 t< pxxxX 2 + 7 G$^xxX + 3G^x)^ — 6 H 2 (5 G^^xX + G$^ 
[2i7 3 (2 G^.^xxX 2 — 9 G^^ — 7 G^^xX) — 12 H(G^.^xX + G^^)H](j) 

6H 2 X (3G 5 ,4>x + 2 Gs^xxX) H + 6 H 2 X (3 H 2 G^^x — £5,000 + 2 H 2 G^^xxX — 2 G^maxX) 



2G 



5,4>4>xx 



X 2 )]i 



'5,<p<p<pX^ 



(35) 



In deriving the above we used the background equations. The expression of the coefficients of the equations (|27)|) - ([2"!)1) . 
written in terms of the variables 9, S, etc.. becomes compact, though many of the coefficients include terms <j) n 
(n > 0) in the denominators. However, there are no divergences at cf> = 0. In fact, whenever the term <p n appears in 
the denominator, the numerator of the same coefficient compensates with the term 4> n . For instance, this property 
can be seen in the expression of M 2 . From the expressions of the coefficients given above, it is clear that not all these 
coefficients are independent. For example, later on, we will find it convenient to use the following relations ^3 = Bs, 
D 7 = B 7 , and D w = A e . 

The mass of the field 5cj) is related with the term M 2 defined in Eq. (f55|) . In fact, for a canonical scalar field 
described by the Lagrangian K = X — V(cf>) with Gi = (i = 3,4,5), we have that M 2 = —K^ = V^. In 
viable dark energy models based on f( R) g ravity and Brans-Dicke theory with a field potential, the term —K^ 
is the dominant contribution to M 2 [43l. |44| . The term —K^ comes from the time-derivative of /j,, such that the 
contribution from the term 3H(fi + v) is usually unimportant relative to /i. 

The equations of motion (|2l^ - ([3"Tj) are not independent. In fact, we find the identity 



k 2 

3HEq, — HEi& — 



p m (E v + 3HE V ) + p m E s = . 



(36) 



This relation can be used in two ways: 1) to check the consistency of the equations themselves; 2) to get some 
equations of motion, which would be missing when some gauge is used from the beginning. For example, in the 
Newtonian gauge (x = 0), the equation -E x | x=0 = 0, cannot be derived directly. However, it is still possible to obtain 
it by using Eq. ([55]). 

We note that the following combination of the perturbation equations is useful: 



k2 E 



3 I EL 



3HE X - = , 



(37) 



which is written explicitly as 



E~, 



5 fi $ + B 7 



B 8 * + B 10 x + B u x = . 



(38) 



This equation corresponds to the traceless part of the gravitational field equations. 



(i 



In order to study the evolution of matter perturbations we introduce the gauge-invariant density contrast 

S m = 5 + 3Hv . (39) 

From Eqs. (J30) and (J3T|) it follows that 

k 2 /•• -\ 

6 m + 2HS m + -(*- x )=3(l + 2Hl), (40) 

where / = Hv — <J>. 

IV. EFFECTIVE GRAVITATIONAL COUPLINGS UNDER THE QUASI-STATIC APPROXIMATION 

ON SUB-HORIZON SCALES 

When we discuss the evolution of matter perturbations relevant to large-scale structure and weak lensing, we are 
primarily interested in the modes deep inside the Hubble radius (k 2 ja 2 ^> H 2 ). We shall use the quasi-static approx- 
imation on sub-horizon scales, 2 under which the dominant contributions in the perturbation equations correspond to 
those including k 2 /a 2 and 6 0, H3, 51 . There 

are two different classes of dark energy models, depending on the mass 

M of a scalar degree of freedom. 

The first one corresponds to the case in which the mass M becomes large in the early cosmological epoch. The viable 
dark energy models constructed in the framework of f(R) gravity [43], 53, 51] and Brans-Dicke theories [44| belong to 
this class. Since the effect of the field mass cannot be neglected in such cases, we need to take into account the term 
M 2 to discuss the evolution of perturbations. This induces the oscillation of the field perturbation d<j>, but as long 
as this oscillating mode is initially suppressed relative to the matter-induced mode, the quasi-static approximation 
can reproduce numerically integrated solutions with high accuracy (43l |49| . Under the quasi-static approximation 
we neglect the time-derivatives of 8<f), which corresponds to the approximation under which the oscillating mode is 
unimportant relative to the matter- induced mode. 

Another class corresponds the case in which the field does not have a massive potential, e.g., Galileon gravity. In 
such cases the numerical simulations in Refs. (35], |38| also show that the quasi-static approximation is sufficiently 
accurate for the modes deep inside the Hubble radius. 

We expect that the quasi-static approximation on sub-horizon scales should be trustable for our general theories as 
well, provided that the matter-induced mode dominates over the oscillating mode. 

Let us choose the Newtonian gauge in which % = 0. Under the quasi-static approximation on sub- horizon scales we 
find that Eqs. (|3"5]h (gjj) , and (J2BJ) can be rewritten as 

B 6 <S> + B 7 5(j) + = 0, (41) 

k 2 ( k 2 \ k 2 

B 7 —<5> + D 9 — - M 2 Scb + A 6 —y ~ , (42) 
a z \ a z J a z 

B 8 ^<f> + A 6 ^r5<t> - p m 6 ~ 0, (43) 



where we used the relations A% = Bs, D 7 = B-j, and D\q — Aq already mentioned in Sec. Mil The above three 
equations correspond to the traceless part of the gravitational field equations, the scalar-field equation of motion, and 
the (OO)-component of the gravitational field equations, respectively. Note that Eq. (|4ip also follows from Eq. (|2T[) in 
the same approximation scheme. 

We can solve Eqs. (|4T|) and (|4"2l for $ and 5(f) in terms of 'J, and then substitute these expressions into Eq. (|4"3"1) . 
This gives the following Poisson equation 

k 2 

— * ~ -4nG eS p m S . (44) 
a z 



2 Strictly speaking, the typical scale here should be given by the sound horizon rather than the Hubble horizon because the propagation 
speed of the scalar mode, c s , differs from unity in general |2ll. l46l|. One needs to be careful for the use of the quasi-static approximation 
in models with c s < 1, as the range of the validity of the approximation may be quite limited. 
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Here the effective gravitational coupling G G ff is given by 

2M 2 l [(B 6 D 9 -B 2 )(k/a) 2 -B 6 M 2 } 

G e ff = o G (45) 

(A 2 B 6 + B%Dq - 2A e B 7 B s ) (k/af - BfM 2 

m p\ { (® + ffe ) + (e?T - eg T /e) 2 + j- T [iM 2 a 2 /t 2 + (f + p)/2] 

'G, (46) 



6 2 T s + G 2 T [XM 2 a 2 /k 2 + {£ + V)/2] 



2 



where G is the bare gravitational constant related with the reduced Planck mass M p \ via the relation 8irG = M pl 
and 

*-§l(S*)-*- («» 

In order to avoid the Laplacian instability of scalar perturbations we require that Tg > [2ll , l4ll . [4^ | . While G e ff 
is written in a compact expression in Eq. (|46|) . it is often convenient to use the form (|45|) for a given Lagrangian. 
In Appendix A we present the explicit forms of the coefficients B§, By, Bg, and Dg, which is useful for the 
computation of Eq. (j^5l) . 

Under the quasi-static approximation on sub-horizon scales Eq. (|40p gives 

k 2 

S rn + 2H6 m + — * ~ , (48) 

On using Eq. (gj) and 5 m ~ 5 (which are valid for k 2 /a 2 3> -ff 2 ), it follows that 

6 m + 2H6 m - AnG cS p m 8 m ~ . (49) 

This can be written as 

C + ( 2 + §) S 'm l^mS m * , (50) 

where Q, m = p m j (3M 2 \H 2 ), and a prime represents a derivative with respect to N = In a. 

We define the anisotropic parameter r] to characterize the difference between the two gravitational potentials: 

»7 = -$/#. (51) 

Under the quasi-static approximation on sub-horizon scales this reduces to 



(B 8 D 9 - A 6 B 7 )(k/a) 2 - B 8 M 2 

(B 6 D 9 - B 2 )(k/a) 2 - B 6 M 2 
g T + {HF T - Q T )Q - HQt{Gt + HQt) - Qt [(£ + V)/2 + XM 2 a 2 /k 2 



(53) 



j- T (e + ire) - (e T + #£t) 2 - t t [{£ + p)/2 + xM 2 a 2 /k 2 ] 

We also introduce the effective gravitational potential 

$ ofr = (tf - *)/2 , (54) 



which is associated with the deviation of the light rays in CMB and weak lensing observations [50j. From Eqs. j44] 
and (1511) we have 



<I>, 



a n 1 ±V( a \ 2 x 3 G cS 1 + t] f aH\ 2 
4^G cff ^-(-j Pm 6^-- — —\ — \ n m 5 m . (55) 



Let us consider k-essence in the framework of General Relativity (GR), which corresponds to the Lagrangian C = 
K((j),X) + {M 2 l /2)R [i.e. G 4 = M^/2, G 3 = = G 5 ]. In this case one has T T = Qt = M 2 t , 6 = #M 2 U and 
■7-s = —M^H/H, which gives G e g = G and rj = 1 from Eqs. (|46|) and (f53|) . During the matter-dominated epoch 
(if' / H ~ — 3/2 and f2 TO ~ 1) there is a growing mode solution J m cx a to Eq. (|5T))) . For this solution $ ff = constant 
from Eq. ()55p . In modified gravitational theories G G ff and rj are in general different from G and 1 respectively, so that 
the evolution of 8 m and $ ff is subject to change compared to GR. 
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V. APPLICATION TO SPECIFIC THEORIES 

In this section we apply our formulas of G e ff and r\ derived in Sec. lIVI to a number of modified gravitational theories. 

A. f(R) theories 

The Lagrangian of f(R) theories corresponds to C = (M pl /2) f(R), where / is an arbitrary function in terms of the 
Ricci scalar R. This is equivalent to the Lagrangian JT]) by choosing the following functions [5l| 

M nl 1 

K = f(Rf.R-f), G 3 =0 = G 5 , G 4 = -M pl f = M p i/,«, (56) 

where is a scalar degree of freedom having the dimension of mass. Since G4 = <pM p \/2, B-j — 2Aq — 2M p \, 
B & = Bg = 2M p y<j>, and D 9 = in this case, Eqs. flU} and ([HJ read 



_ M pl 4 + 2(0/M pl )(Ma/fc) 2 _ 1 + (cf>/M pl )(Ma/k) 

° ff 3 + 2(0/M pl )(M a /fc) 2 ' 1 2 + (<f>/M p i)(Ma/k) 2 ' 

From Eq. (|35p the mass squared of the scalar degree of freedom is given by 

1 



M 2 = -K M = — — . (58) 

4J,RR 

Substituting this relation and <fi = M p \J.r into Eq. (|57p . it follows that 

_ G l + 4(/, flfl //, fl )(fc/a) 2 _ l + 2(/,W/,«)(fc/") 2 r , ql 

Cff /,* 1 + 3(/,W/,fl)(fc/a) 2 ' l + 4(W/,fl)(W 

which agree with those derived in Ref . [3l| . 

The viable dark energy models based on f(R) theories were constructed to have a large mass M in the deep matter- 
dominated epoch [4^, [43, [4^| , i.e. J,rr 3> 1 and Jr ~ 1 for R ^> Hg, where Ho is the Hubble parameter today. In 
the regime (f,RR/f,R.)(k/a) 2 <C 1 (or M 2 f^R ^> k 2 /a 2 ) one has G e s — G and 77 ~ 1, so that the evolution of density 
perturbations is similar to that in GR. At late times the mass term M gets smaller with the growth of J.rr. Since 
G e fi — 4G/{Zf,n) and 77 ~ 1/2 for {f ,rr/ f ,R){k / a) 2 1, the growth rate of matter perturbations is larger than that 
in the ACDM model, e.g., S m oc t(v / 33~ 1 )/ 6 during the matter-dominated epoch [H, |H|. 

Substituting Eq. (|59|) into Eq. (|55f we obtain $ e ff ~ —(3/2)(aH/k) 2 £l m S m , where Cl m = £l m /f,R- This means 
that the anisotropic parameter 77 between the two gravitational potentials practically compensates the modification 
induced by the gravitational coupling G c g, i.e. G e g (1 + r))/2 = Gj f r. During the matter dominance (0 TO ~ 1 and 
a oc r; 2 / 3 ) the evolution of the effective gravitational potential on sub-horizon scales is given by $ e ff oc ^(v / 33-5)/6 "4<| 

B. Brans-Dicke theories 

Brans-Dicke theories [l(| with the field potential V{4>) correspond to the choice 

K= Mpl ^ X -V(4>), G 3 =0 = G 5 , G 4 = ±M pl <t>, (60) 

where wbd is the Brans-Dicke parameter (which is constant). Compared to original Brans-Dicke theories we have 
introduced the reduced Planck mass M p \ in K and G4 such that the field <fi has a dimension of mass. The difference 
from f(R) theories appears for the kinetic term M p \wby}X / <p, in which case Dg = —M p \uibd/<P- From Eqs. (|45[) and 
([5*2"]) it follows that 

r _ Mgl 4 + 2uj bd + 2(<j ) /M pl )(Ma/k) 2 _ I + a; B p + (<j)/M pi )(Ma/k) 2 

cff 3 + 2 WB D + 2(0/M pl )(Ma/fc) 2 ' ?? 2 + cjbd + (<t>/M p \)(Ma/k) 2 ' 1 J 

where 



M 2 = V*. + 



(0 + 3i^)] . (62) 
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The results (|57|l in f(R) theories can be recovered by setting wbd = in Eq. (f6"Tj) . It is convenient to express 
M 2 solely in terms of the potential. This can be done by using the scalar field equation of motion, and we obtain 
M 2 ~ + V^/(f>, where we have neglected 0(M P \H 2 /4>) terms. 

In the limit where wbd — > oo or M 2 — ► oo (with <fi ~ M p \) we recover the General Relativistic behavior: G e ff — G and 
■q ~ 1. In the limit that M 2 -> we have G off ~ (A/ pl /</0(4 + 2w BD )/(3 + 2w BD )G and r) ~ (I+w BD )/(2+w BD ) [H,|3l|. 
The effective gravitational coupling in the latter case also agrees with the one corresponding to the gravitational force 
between two test particles [52| . 

C. Kinetic gravity braidings 

Let us consider the kinetic gravity braidings [53[ described by the Lagrangian 

1 

2 



K = K(<f>,X), G 3 = G 3 (0,X), G 4 = ^M p 2 l5 G 5 = 0. (63) 



Since A 6 = — 2XG 3i x, B 6 = B s = 2M 2 V B 7 = in this case, it follows that 



r M 2 -D 9 (k/a) 2 

^ cff - M 2 -(D g + 2X 2 Gl x /M 2 1 )(k/a) 2 °' ?? ^ 



where 



# 9 = -A',* - 2 (G 3) x + ^G 3 ,xx) 4> - 4ffG 3 ,x^ + 2G 3 , - 2XG 3 , 0X . (65) 

In the limit that M 2 — > oo we have G e g — ?> G, so that the General Relativistic behavior is recovered. 

Let us consider the theories in which both K and G 3 depend only on X, i.e. K = K(X) and G 3 = G 3 (X). Since 
7\/ 2 = 0, G 3 .0 = 0, and G 3i ^x = in such theories, the effective gravitational coupling is given by 



G c ff = G < I + = ^ r- ^ . (66) 

\ + 2(0 + 2i70)G 3 ,x + G 3 , XJf 2 </)] - G 2 ^ 4 J 



This result agrees with that derived in Ref. [37[ in which the authors studied the evolution of perturbations for the 
functions K = —X and G 3 oc X n (which corresponds to the Dvali and Turner model (23|). 

One can also extend the analysis to the case where G4 is a function of <f>, i.e. £ = Gi{4>)R + K((j>, X) — G 3 (<$,X)D<A. 
The perturbation equations for the theories with Gs((/>,X) = £(4>)X were derived in Ref. [36[ (see also Refs. [34l l35| 
for specific choices of the functions G±(4>) and K(cj),X j). 



D. Covariant Galileon 

The covariant Galileon without the field potential corresponds to [54| 

K = -c 2 X, G 3 = ^X, G^hl^-^X 2 , G 5 = ^X 2 , (67) 

where c, (i = 2,3,4,5) are dimensionless constants and M is the constant having the dimension of mass. For the 
choice (|57|) we confirmed that all the coefficients in Eqs. (|32p - (j33)) are equivalent to those given in Ref. [Hj]. Hence 
Eqs. (|45[) and (|52p reproduce the effective gravitational coupling G e s and the anisotropic parameter 77 derived in |38l| . 

For the covariant Galileon there exists a stable de Sitter solution where X = constant [55[ (see also Refs. (56ll57| for 
related works). Since G e ff is larger than G before the solution reaches the de Sitter attractor, the growth rate of matter 
perturbations is larger than that in the ACDM model. In addition the variation of the effective gravitational potential 
^eff can be more significant than that in f(R) gravity, because 77 can be larger than I in the early cosmological epoch 

m 

E. Field derivative couplings with the Einstein tensor 

The dark energy model of Gubitosi and Linder (26| corresponds to 

K = X, G 4 = ±M 2 ,, G 5 = -A^, (68) 



10 



where A is a dimensionless constant (see also Refs. |25l . |58|). The Lagrangian in Ref. [26j involves the term 
(A/Mp^G^V^V 1 ^, but this is equivalent to -A^/M^G^V^V 1 ^ after the integration by parts. 

Since A 6 = -AXH^M^, B 6 = 2M^ + AXX/M^, B 7 = -4A(0 + H<j>)/M^, B 8 = 2M^ - AXX/M^, and D 9 = 

-1 - 6XH 2 /M^ - iXH/M^, it follows that 

Goft = fl + |jG, (69) 

2X[M 2 l <f> 2 + 2\{2H<j>4> + 24? + ZH 2 4> 2 + 2H^ 2 )] 
V ~ 1 M° + XM^MI^H 2 + 4M^H + 4> 2 ) + 2\ 2 {7H 2 4> 2 + 2H4> 2 + 40 2 + 8H<fxf) ' ^ 

where 

£1 = A[3AfpV^ 2 + AA/p\(18A/pVff 2 2 + 1^^^ ^ 
6 = Mpi + 2\M% l (3M 2 l H 2 + 2M^H - <j> 2 ) + X 2 M^ 4 + WM^H^ - AM^H 2 ^ 2 - 8M^H^ 2 ) 

+2\ 3 <fi 3 (2H<j)-8H4>-9H 2 <i)) . (72) 

For A ^ one has G e s ^ G and rj ^ 1, so that the evolution of perturbations is different from that in GR. Note that 
the results ([55)) and ([70|) are derived for the first time in this paper. 



VI. CONCLUSIONS 



In this paper we have derived the full equations of scalar density perturbations for the perturbed metric (|2ip. The 
Newtonian gauge corresponds to the choice x = 0, which fixes the temporal part of the gauge-transformation vector. 
Since the different gauge choices (such as S(j) = 0) are possible, our linear perturbation equations can be applied to 
other gauges as well. 

For the perturbations deep inside the Hubble radius the quasi-static approximation employed in Sec. HVl is accurate 
for an effectively massless scalar field. In fact this is the case for kinetic gravity braidings and covariant Galileons. In 
the dark energy models based on f(R) theories and Brans-Dicke theories, the mass of the scalar field degree of freedom 
needs to be large in the region of high density for consistency with local gravity constraints. In order to accommodate 
such cases we have taken into account the effective mass M in estimating the effective gravitational coupling G e g. 
For the quasi-static approximation to work it is necessary that the time-derivatives of the field perturbation 8(j> 
are neglected relative to the terms including c 2 {k 2 / l a 2 )5(fi, where c s is the scalar propagation speed whose explicit 
expression is given in Ref. (46| . This implies that the perturbation field cannot be fast oscillating. In other words, 
provided that the oscillating mode of the field perturbation is suppressed relative to the matter-induced mode, the 
quasi-static approximation on sub-horizon scales can be trustable. 

In order to estimate the growth rate of perturbations relevant to large-scale structure and weak lensing, it is 
sufficient to use the approximate results of (|45|) and (|55|) with the anisotropic parameter 77 given by Eq. (|52")h We 
applied our formulas to a number of modified gravitational models of dark energy and found that they nicely reproduce 
the previously known results. There are some new models-such as the field derivative couplings with the Einstein 
tensor-in which the evolution of perturbations deserves for further detailed investigation. 

For the large-scale perturbations associated with the integrated-Sachs- Wolfe effect in the CMB anisotropics the 
sub-horizon approximation in Sec. II VI is no longer valid. Instead we need to integrate the perturbation equations 
(|26]) - (|3ip numerically, along the lines of Ref. [38|. It will be of interest how the joint data analysis of CMB combined 
with the observations of large-scale structure and weak lensing place constraints on each modified gravitational model 
of dark energy accommodated by the general action ©. 
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Appendix A: Coefficients relevant to sub-horizon perturbations 

Here we write the explicit forms of the coefficients A$, B$, B7, B$, and Dg: 

A 6 = -2XG 3 ,x-4H(G A ,x+2XG^xx)<P + 2G^ + 4XG^ x 

+AH (G 5 , + XG^x) 4> - 2H 2 X (3G 5 , X + 2XG 5 ,xx) , (Al) 
B, = 4[G 4 -A(0G 5 ,x + G 5 , )], (A2) 
B 7 = -4G A , x H<j>-4(G iiX + 2XG± } xx)4> + 4G4 1 <i>-%XG^x 

+4(G 5 , + AG 5 , 0X )0 - 4H[{G 5jX + XG 5 , X x)4> - G^ + XG^ X ]4> + 4A[G 5 , 00 - (H 2 + H)G B ,x], (A3) 
B s = A[G i -2XG 4>x -X(Hj>G 5t x-G 5t t)], (A4) 
D 9 = -K x -2(G 3 ,x + XG 3 ,xx)4>-4HG 3 ,xj> + 2G 3 ,4 > -2XG 3 ^x 

+ [-4 H(3 G^xx + 2XG i ,xxx)4> + 4# (3G 4 , x - 2XG^ X xM + (6 G^ x + *XG^xx)4> 

— 2QH 2 XG^xx + ^XGi^x — ^H{G^x + ^XG^^xx) — &H 2 G^x 

+{4ff(2G 5 , x + XG 5 ,txx)4> - *H[(H 2 + H)(G 5 , X + XG 5}X x) ~ XG 5 ,^x]}4> - 4H 2 X 2 G 5 ^ X x 
-2H 2 (G 5 , X + 5XG 5 , XX + 2X 2 G 5 , X xx)4> + 2(3i/ 2 + 2H)G^ + AHXG 5 ^ X + 10H 2 XG 5 ^ X - (A5) 
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